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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM 1336 


DEVELOPMENT OF A LAMINAR BOUNDARY LAYER 
BEHIND A SUCTION PoINT* 


By W. Wuest 
1. INTRODUCTION 


Boundary-layer suction originally was applied to reduce the 
boundary-layer thickness and therewith the inclination to flow sepa- 
ration; however, since the properties of bodles with small drag have 
been improved more and more, attention was drawn to an increased extent 
to the reduction of surface friction. One now strived toward keeping 
the boundary layer laminar as long as possible, thus to defer the tran- 
sition point to turbulence as far as possible. DBoundary-layer suction 
was recognized to have a favorable effect in this sense, and therewith 
the velocity distribution in a laminar boundary layer behind a suction 
point acquired heightened interest. The stability of a laminar velocity 
profile is very severely affected by the shapé of this profile. = 2 © 

In a considerable number of theoretical reports (reference 1) the 
case of continuous suction was treated for reasdéns of mathematical 
simplicity; permeability of the wall surface was assumed. In further a 
reports, the stability of laminar boundary-layer profiles in case of 
continuous suction was treated and a considerable rise in the stability 
limit was determined; however, a technical realization of such perme- i = 
able walls with sufficiently smooth surface and adequate material 
strength characteristics is difficult. For structural reasons, it is 
simpler to arrange single-suction slots. In addition to the suction 
effect proper, there appears here the sink effect first discussed in 
detail by L. Prandtl and O. Schrenk (reference 2) and recently treated © 
by Pfenniger (reference 3) in an instructive experimental investigation. 

Below, the pressure variation along the wall as well as, in partic- 
ular, the sink effect are disregarded. Figure 1 shows the practical 
realization of such a case. We assume that on a flat plate A, a laminar 
boundary layer ("Blasius boundary layer") develops at constant pressure. 
We assume a second plate B arranged beginning from a certain point x, 


at the distance y, parallel to the first plate so that a suction slot 





"Entwicklung einer laminaren Grenzschicht hinter einer 
Absaugestelle,” Ingenieur Archiv, Vol. 17, 1949, pp. 199-206. | Br 
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is formed between the two plates. The magnitude of the power require- 
ment for suction is assumed to be precisely such that merely the part 
of the boundary layer situated between the two plates is removed. ‘Thus, 
there begins above the plate B a new laminar boundary layer which is. 
distinguished from the Blasius boundary layér by another initial condi- 
tion. The new boundary layer forms at its start the outer part of a 
Blasius boundary layer. OO 


2. BOUNDARY-LAYER EQUATION AND ASYMPTOTIC BEHAVIOR 


By introduction of the stream function and the total pressure, the 
boundary-layer equation may be transformed by the well-known method 
(reference 4) into 


> En N 
og Og u 
Ox age oe 


where g =p + Şu u? and = u. We limit ourselves to the case thet 
the flow takes place outside of the boundary “layer at the velocity 
u] = const., thus to the flat plate and put furthermore 


g =-6 u 1 - a(x,¥)) + Conste i (2) 


or, respectively 
u = u, ya Rn | sh (3) 


This statement has been chosen so that for large y-valués, q assumes 
the value 1. Equation (1) is thereby transformed into v 


2. . 
Od _ vu, ya 24 | (4) 


From the definition of the stream function and from equation (3), one 
further obtains with n = WA /va, x ' 
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du _ 1 fee) Sa (6) 


Ka ne 
y=x el, ia (7) 


In order to investigate the asymptotic behavior of the differential 
equation (3), we put for large values of y 


q=l-q, with q«l (8) 


In first approximation, one then obtains 





= yu] u | (9) 


This differential equation, however, is mathematically identical with 
the differential equation of a nonsteady flow independent of x which 
has been treated before (reference 5); the time t is now replaced by 
the stipulated space coordinate x. It also corresponds to the well- 
known heat-conduction equation. The general solution is therefore given 
by s _ 


x) = Ł HEN N 6 Vig 
hiki a ny si (Pac a: ) (rac - “xy u 


0) 


(0,28 of —__¥ — \üxt (10) 
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Therein 


zZ D - 
Velo 


is the knom error integral. W. Tollmien (reference 6) has investigated 
this solution for two special cases where the first integral disappears. 
For the boundary layer with suction, however, this will no longer be the 
case. 


3. BLASIUS BOUNDARY LAYER 


Although we presupposed that the velocity u] at the edge of the 


to be treated here nevertheless differs from the flow on a simple flat 
plate ("Blasius boundary layer") by the fact that other_initial condi- 
tions exist; rather, the Blasius boundary layer is contained as special . 
solution among the suction boundary layers since there "X, = O, thus 


suction point and beginning of the plate A (fig. l) coincide. Since we 
shall make use of this special solution for the later calculation, we 
shall first consider the Blasius boundary layer. It is distinguished 
by the fact that q may be regarded as dependent merely on a quan- 


tity n = y} vax. One then obtains from equation (4) the following 
differential equation of the Blasius boundary layer | 





Sie aay nn 
i E 5 à k #0 a = (11) 
n on 


The solution may be written in the following form 


3\1/2 3\2/2 3\3/2 
ee) om 
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The constants a, therein have the following values 


a ee a5 = 0.57627 x 107° ag = 3.7282 x 10°“ 
BERN 22.8 or 2 -13 
Bp == 55 a, = 3.8907 x 10 3 ajg = 2.2383 x 10 
= -T u cde -9 = , 10l} 
Sy. ye E a7 = -1.3986 x IEO a, 7 = -0.3104 x 10 
3 (990 x 15) 2 ' = 
a, = 1.60333 x 107° ag = -3.9135 x 10-11 ap = -1.081 x 10 


Due to the boundary condition at thg wall, one integration constant is | 
zero. The second integration constaft is determined from the asymptotic 
behavior for large values of n. Because of qg(¥,0) = O the first 


integral in equation (10) is elimingted. The second integral, however, 
yields by partial integration, with consideration of the asymptotic 
behavior of the error integral, just as in the case treated before by 
W. Wuest the solution Š . 


i 
wiz]. - of (13) 
-E í i Ken) | 


The constants ß in equation (12) and y in equation (13) are deter- 
mined by the fact that for large n values q and daq/dn according 

to equation (12) and equation (13) agree yith each other. The recalcu- 
lation of the two constants yielded the following values = =~ — 


B pO. 6642 `. -y 7# 0.828 
> > 
y 
For comparison, L. Prandtl (reference 7) gives the P6flowing values 
calculated by Blasius ang TREE oe which read, converted to the above 
designations oe 


2 x 0.332 = 0.664 


Ww 
1i 


av 0.231 


~“ 


.0.819 z a 


`g 
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F. Riegels and J. A. Zaat give ina new report, (reference 8) for y 
the following value 


y = 0,342 yen = 0.857 


The function q with first and second derivative has been tabulated 
and plotted in numerical table 1 and figure ol, 


4, ASYMPTOTIC BEHAVIOR OF THE SUCTION BOUNDARY LAYER 


For calculation of the asymptotic behavior of the suction boundary 
layer, we divide the function q, defined by equation (8) into two 
parts 


ty = Sy * Tye u 


The first part is to be selected so that it satisfies the initial condi- 
tion at the suction point x = Xg; this is done. by extending the asymp- 


totic solution of the Blasius boundary layer to x > Xọ as well. From 
equation (13) one then obtains 


© Ati 
dyi = 7t - i 
hvu; x 


lhe numerical table has been calculated with the values ß = 0.664 
and y = 0.819. 


“4 


NACA M 1336 T 


Numerical Table 1. Blasius Boundary Layer 


n a(n) a(n) at? (n) yaln) 
0. 6640 
.1 . 6555 
.2 . 6427 
3 . 6206 
ot . 5981 
«5 25734 
.6 ‚5471 
“7 «5195 
8 4912 
.9 . 4625 
1.0 4337 
1.1 . 4062 
1.2 . 3766 
153 . 3487 
1.4 . 3217 
19 . 2959 
16 . 2705 
LT . 2466 
1.8 .2238 
1.9 .2022 
2.0 . 1820 
2.5 .1013 
3.0, 0509 
35° 0217 
1,0 .0085 
4,5 0029 
5.0 0009 





Therein y= 0 forms the new wall streamline and Yọ the suction 
Quantity. The second part d;o then must be chosen so that the bound- 
ary condition q,=4,(0,x) is satisfied. If the asymptotic rela- 
tion q~ l =- qs would rigorously apply in the entire domain of the 
boundary layer, there would have.to be at the wall q,(0,x) = 1, because 


of q = 0; however, the asymptotic solution deviates from the rigorous 
solution if it is continued up to the wall. Therefore qd; (0, x) is an 


unknown function regarding which we merely make the assumption that it 
does not become infinite. As initial condition for the part q,n one 


further has q.,0( Vs XQ) =O since q.1( ¥ Xo) 
initial condition 


\ 
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already satisfies the 


which insures the connection with the Blasius ‘solution. The contri- 
bution qo to the solution alb must obey the differential equa- 


equation (10)) the first integral is elimi- 


tion (9). 


In the solution ( 
nated, because of a,o(¥, 


has to put 


4,o(0,x) = 4,(0,x) - q,(0,x) = a (0,x) - 7/2 - Gs) 
| hvu, x 


so that the asymptotic solution reads 





ee ah 
vu x Xo 


By partial integration one obtains with considerati 
behavior of the error integral 


Y 
dy 7jl - ( 


(9 %0) 


so that one finally obtains as the as 


boundary 


+ % 





Vw, x 


Because of the comnection with the Bla 
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Xo) = 0, whereas in the second integral one _ 


on of the asymptotic 


(by W. Wuest, elsewhere) 


Y 


sius solution, however, 





en) 


= y, if the asymptotic solution is continued up to the wall, 


layer 


L 


ymptotic. solution for the suction 
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Instead of the error integrals © one may for large values of y 
again go back to the Blasius solution if one takes the asymptotic 
behavior of the latter according to equation (8) and equation (13) into 
consideration 


pers ah Gow AB E E De De (15) 
Vu, x V4Vu xo yvu (x - Xo) 


In this formula ap represents the Blasius solution. The last form of 


the solution proves to be particularly expedient for the further con- 
siderations. 


5. APPROXIMATE SOLUTION FOR THE SUCTION BOUNDARY LAYER 


It suggests itself to generalize the asymptotic solution which is 
valid for large values of wy in the following manner 


sal oeiia iea cee ARREHN (16) 
A) EEN 


Due to q=0 for y=0 and because of equation (15) the func- 
tion F(y,x) must satisfy the following conditions 


| y y 
B ER 
vux Vu} Xo 


It was hoped at first that one could choose for F, as in the nonsteady 
analogue by W. Wuest, elsewhere correspondingly an exponential func- 
tion as the simplest formulation; besides equation (17) thé disappear- 
ance of the second derivative of q at the wall would be added as a 
further condition; however, it was shown that such a formulation does 
not meet with success and even, in a certain domain, does not yield any 
solution at all. 





(17) 


-.: 


10 NACA TM 1336 


For the further calculation we introduce the following simplified 
notation u 





y+ Y Y | X = | 
BEE =f l = ni! O = No n = No + ——— n! (18) 
vu, x vu, (x xo) yvu]x 


so that the solution (equation .(16)) reads 


a= q,(n) - F| (2 : a(n") (19) 


According to a suggestion by A. Betz, we equate as first approxi~ 
mation 4 the function F to the value dependent only on x 





Yo 
Fo(x) = F(0,x) = af =) = (No) (20) 


at the wall. Thus the first approximation reads 


pn) - Fol (2 - 0) | (21) 


This formulation does not fulfill the boundary-layer equation (4) exactly. 


In particular, the second derivative of qy “at the wall does not disap- 


pear; however, the dependency on the second derivative of the stability 
of the velocity profile is of a very sensitive nature so that one has to 
look for a more accurate solution. By substitution of the approximate 
solution {equation (21)) into the boundary-layer equation (4), one 
obtains 


tta 0 ae t = 3 32e 
D = RE- 1% i} = 2( ~ %) fay ‘ro 


_ 


Hence there results oe, font as the error of this First approxi- 


mation. By subtraction of the exact solution in which F stands 
for Fọ and q for 4, while Gy disappears, one then obtains 


en 


EL 
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ao leer [2 - (ni fep 2 
-To n (22) 
Sn! 2 ) da ‘<7 "a2 
where 
is an unknown function. The quantity G5! ' disappears for n’ = 0 
and yt = % By integration of equation (22) one obtains 
(23) 


(Œ - Fo) [1 - ana] ee | 


o be determined graphically or numerically by repeated 


Therein 67 is t 
quadrature 
| (2h) 


No 3° 
u -f T a = 
nn 


Ep by substituting in the 


the asymptotic 


Thereby one obtains 


mine the asymptotic behavior of 


J! for \/a, and Va 


i 
2 dr 


One may dete 
above definition of 60 


values ya ~ 2-3 4a and q~i- 
eot! ~ (x - xo) S(r- 1) 
ox 

ions (9), (19); (21), and repeated inte- 


Hence follows with use of equat 


gration with respect to n? = yva (x - Xo) 


Ep ~ (Fr - Fo)|} - (| (25) 
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As before, © denotes the error integral. Generally we visualize «o 
as represented in the following manner -= 


+>. a(x) - (+) | (26) 


HL 


By way of approximation we limit ourselves to the first two terms, 
with a, = 7(F,- Fo) and a determined by the fact that q must 


disappear at the wall. We determine accordingly the function F 
approximately to be 


F = Fo + een) pe + 1( Fo - Fo) ul (3) + a, ea} 
(27) 


a, = ~6, (0,x) - Y(Fx ~ Fo) (28) 


Calculation example.- Vo = 0.125 was selected as numerical 
example; F was calculated for the values x/xp = 1.234, 1.562, 4.34, 
and 9.78 and plotted in figure 3. For x/Xp = 1.562 the error was 
determined by substitution of the approximated solution into the 
boundary-layer equation, and compared with the first approximation 
according to equation (21). Compared to the first approximation, a con- 


siderable improvement results particularly in the region near the wall 
(fig. 4). In figure 5 the results are converted to the velocity pro- 


file, in figure 6 the second derivative is represented. As a supplement, 


the connection between the degree of suction and the suction quantity of 
the magnitude 1,* = Wo ES will be suppletiented. By the degree of 


suction @ we here understand 


osla  - 7 (29) 


5,* being the displacement thickness immediately ahead of the suction 
point and 6 5* immediately behind it. Therewith @ is given by 


a 


è 
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(30) 


The resulting values are tabulated in table 2 and plotted in figure 7. 


Table 2. Degree of Suction 


ox [or Tor [os [os [oo [oo 


The suction quantity Vo is, furthermore, given by the following 
relation 

















ua Ur VE (31) 
6. SUMMARY 


ö ~ 3 
The rable gpg d laminar boundary layer behind a suction point 
is investigated if by¥f the suction merely the part of the boundary layer 
near the wall is "cut off", without the slot exerting a sink effect. 

As basis of the calculation, we used the boundary-layer equation in the 
form indicated by Prandti-Mises which is.closely related to the heat 
conduction equation or, respectively, to the differential equation of 
the nonsteady flow(which is independent of the coordinate x along the 
wall. With consideration ea E behavior of the solution, 
an approximate solution is developed which is similar in structure to 
the solution of the nonsteady analogue which has been treated in an 
earlier report by W. Wuest,.elsewhere. | 


Translated by Mary L, Mehler 
National Advisory Committee 
Por Aeronautics 


Im 
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Figure 2.- The function q(n) of the Blasius boundary eyen with first and 
second derivative. - 
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Figure 3.- Auxiliary function F(n‘,x) for calculation of the suction 
boundary layer for Vo/yvuız = 0.126. 
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Figure 4,- Error of the first and second approximation for -x/%, = 1.562. 
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Figure 5.- Velocity profiles of the suction boundary layer for 


vof vVuXo = 0.125 and various distances from the suction 
point. 





O 0,050 0.100 0,150 


Figure 6.- Second derivative of the velocity profiles of the suction 
boundary layer for to / vux = 0.125. 
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Degree of suction ® 


— 





Q 02 04 06 08 10 12 14 


Figure 7.- Degree of suction (ratio of the cross -hatched and the total 
shaded area in fig. 1). 
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